A quest for frustration driven distortion in Y2M02O7 
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We investigated the nature of the freezing in the geometrically frustrated Heisenberg spin-glass 
Y2M02O7 by measuring the temperature dependence of the static internal magnetic field distri- 
bution above the spin-glass temperature, T g , using the pSR technique. The evolution of the field 
distribution cannot be explained by changes in the spin susceptibility alone and suggests a lattice 
deformation. This possibility is addressed by numerical simulations of the Heisenberg Hamiltonian 
with magneto-elastic coupling at T > 0. 

PACS numbers: 75.50.Lk, 75.10.Nr 



A geometrically frustrated magnet is one whose sym- 
metry precludes every pairwise interaction in the system 
being satisfied simultaneously. Such systems can have 
"macroscopic" degenerate ground states, which prevents 
the system from settling into a unique state as the tem- 
perature is lowered. In this situation, small perturbations 
to the Hamiltonian can have a dramatic effect, selecting 
one ground state over another. It has been suggested 
that the ground state degeneracy of the Heisenberg anti- 
ferromagnet on the pyrochlore lattice, a lattice of corner- 
sharing tctrahedra, can be lifted by a magneto-elastic 
couplingi*2iiMiiii&. Such a coupling allows the lattice to 
distort in order to relieve the magnetic frustration, con- 
comitantly lowering the total system energy. Since the 
3D corner-sharing tetrahedra lattice is susceptible to dis- 
tortion under arbitrarily small magneto-elastic coupling, 
it is at the centre of an ongoing quest to observe this 
phenomena. Several materials show a lot of promise in 
this respect, including several Cr 3+ spinels and several 
A 2 B 2 7 pyrochloresSiSi 

One promising, but puzzling, candidate is Y2M02O7 
[YMO]. YMO is a narrow-band semiconductor with a 
Curie- Weiss temperature of ~ 200 K and a spin-glass 
transition at 22.5 KS. The glassiness of the system 
implies some kind of disorder in the exchange integral, 
whose origin might be a deformed lattice. Indeed, Booth 
et al£ showed, using the X-ray-absorption fine-structure 
technique, that the Mo tetrahedra are in fact distorted at 
the local level by about 5%. This relatively large amount 
of bond disorder is not seen by the usual diffraction tech- 
niques (X-rays or neutrons), indicating that the aver- 
age bulk structure is almost the perfect oxide pyrochlore 
lattice 8 -. Further evidence of disorder was revealed by 
89 Y NMR™, where regularly spaced peaks implied the ex- 
istence of many non-equivalent 89 Y sites. These features 
were attributed to distortions in the Mo-sublattice. How- 
ever, the NMR data were limited to nitrogen tempera- 
tures due to limitations on the maximum RF power avail- 
able in Helium atmosphere, and the decrease of signal in- 



tensity with decreasing temperature due to the abnormal 
increase of the line width. To gain a better understand- 
ing of the freezing process in YMO, and to determine 
whether a lattice deformation is active close to T g , we 
performed a study similar to the NMR one and with the 
same sample, using muon spin rotation and relaxation 
techniques. In this letter we present these /LtSR results 
and attempt to reconcile them (and the NMR) with the 
existing diffraction data through numerical simulations. 

Transverse [TF] and longitudinal field [LF] /iSR mea- 
surements were performed on the GPS spectrometer at 
the Paul Scherrer Institute, Switzerland. The measure- 
ments were carried out with the muon spin tilted by 50° 
relative to the direction of the applied magnetic field 
of 6 kG, and positron data were accumulated in both 
the forward-backward (longitudinal) and the up-down 
(transverse) directions simultaneously. This allows us to 
determine both transverse and longitudinal muon spin 
relaxations rates. In Fig. ^ we show the LF [panel (a)] 
and the TF data [panels (b) and (c)] at two tempera- 
tures. Both TF and LF relaxation rates increase with 
decreasing temperature; however, at 23.2 K the trans- 
verse asymmetry is zero after 0.6 //sec while the longi- 
tudinal asymmetry is still finite at much longer times. 
This means that close to T g the TF relaxation is larger 
than the LF one and most of the in-plane depolarization 
is from static field inhomogeneities. 

We found that the /iSR LF asymmetry is best de- 
scribed by the root exponential 



A LF {t) = A exp(-{R LF t)i) + B g 



(1) 



where the parameters A arc set by taking into account 
the tilt of the muon spin relative to the longitudinal mag- 
netic field, Rlf is the longitudinal relaxation rate, and 
t is time. Similarly, the TF asymmetry is best fitted by 
a root exponential superimposed on a cosine oscillation 
with a transverse relaxation rate Rtf- 



Atf{i) = A exp(-(R T Ft) 2 )cos(ujt + 
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FIG. 1: (a) Longitudinal field asymmetries for T = 23.2 K 
and T — 45.2 K; (b) and (c) real and imaginary transverse 
field asymmetries for T = 45.2 K and T = 23.2 K respectivly. 
The solid lines are fits to Eqs0and|2] 



Aq, the muon rotational frequency (u> = j^Htf), and 
the background (B g ) are common for all the TF data 
sets. Since the TF relaxation is a result of both static 
field inhomogeneity and dynamically fluctuating fields, 
and the LF relaxation stems from dynamic fluctuations 
only, we can extract the rates of the two processes. The 
static field distribution alone will lead to a polarization 
function of the form 

Pstaudt) = Po exp{-(At)^ 2 )cos{cjt) (3) 

where A = [R^p — R^p] 2 assuming isotropic fluctuations. 

In Fig. |2Jwe plot Rtf, Rlf, and A on a log scale vs. T 
on the upper part of the abscissa. As the temperature is 
lowered the difference between Rtf and Rlf increases. 
In addition A grows exponentially on this semi-log plot 
upon cooling. In the inset of Fig. [5] we present A for 
various fields at T = 22.5 K. One can clearly see that A 
grows with increasing field, and therefore we performed 
the experiment in the highest field available. 

The static relaxation rate A is related to lattice de- 
formation via the muon coupling to its neighboring spins 
and the system's susceptibility. To demonstrate this rela- 
tion let's assume for simplicity that the muon is coupled 
only to one electronic spin S (extension to multiple cou- 
plings is trivial). In this case the magnetic field experi- 
enced by the muon is a sum of the external field Htf and 
the field from the neighboring electron Hj n t = y/^sAS 
where g is the g factor, and ps is the Bohr magneton. 
Here A = A(r) is the coupling between the muon and 
electron spin, which we assume depends on the distance 
between them, although at present the origin of this cou- 
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FIG. 2: A plot of the static (A), longitudinal (Rlf) and 
transverse (Rtf) relaxations rates, on a log scale, versus sus- 
ceptibility and temperature. In the inset we present A versus 
magnetic field at a constant temperature of 22.5 K. 



pling is unknown. In a mean field approximation, S — > 
(S) = %Htf/ 9I^b- Therefore, the muon experiences a 
magnetic field B = (1 + Ax)Htf- Assuming that the 
susceptibility and the couplings are isotropic, the time 
evolution of a muon spin due to static field inhomogene- 
ity is given by: 

Pstauc(t) = J Po cos[ 7ai (1 + A x )H TF t]p(A)dA. (4) 

From Eqs. © and (© we find that p(A) = 
f(5A/{2n\A\})/(2\A\) where the effective width of the 
distribution, 5 A, is given by 6 A — \ A/ {xi^Htf)\ and / 
can be expressed in terms of trigonometric and Fresnel 
functions^. 

If the ratio between A and \ remains constant as the 
temperature changes, the distribution of couplings p{A) 
is temperature independent. Unlike most magnets, this 
is not the case for YMO. In Fig. we also present on the 
lower part of the abscissa the susceptibility measured at 
6 kG with the temperature as an implicit parameter. We 
see that A is not proportional to and not even a linear 
function of \. We thus conclude that the distribution of 
coupling constants broadens extremely fast upon cooling 
below 40 K due to random lattice distortions similar to 
those seen in the NMR work at higher temperatures^. 

The local lattice distortions or strains discussed above 
have been elucidated from magnetic probes (/^SR and 
NMR) in a finite field. At first sight these contradict 
previously published diffraction data where no structural 
deviations from the Fd3m space group were observec&ii. 
More recent neutron powder diffraction results in finite 
applied fields also reveal no unusual thermal dependence 
on the oxide pyrochlore structure^. However, Booth 
and coworkers did find a 5% distribution in Mo-Mo bond 
lengths in zero field. This quandary, the apparent con- 



3 



flict between the data sets, required revisiting the theory 
of magneto-elastic deformation of the pyrochlore lattice. 
We have tried to do this through computer simulations 
where we posed the questions: (I) Can the computer 
find a classical ground state with a lower energy than 
that suggested theoretically? (II) How good are the the- 
oretical approximations? and (III) Since experiments are 
done at finite temperature, what is the effect of thermal 
fluctuations on spin and spatial correlations? 

The Hamiltonian we simulate is the Heisenberg anti- 
ferromagnet with a magneto-elastic term as illustrated 
in the inset of Fig. [31 for a single tetrahedon. The rigid 
bonds between atomic positions of length R are replaced 
by springs whose elastic energy is zero when R = 1 . This 
Hamiltonian is given by 

H = JJ2 Si • S 3 - + J'J^iSi ■ S 3 )5 nj + k - J2(6 nj ) 2 (5) 

i>j i>j i>j 

where fry = |r< — rj\ — 1 is the change in bond length 
between spins i and j, J' is the derivative of the exchange 
J relative to bond length, and k is an elastic constant. 
The minimum energy of TL per spin is given by 



where N is the number of spins 3 ". This energy is achieved 
in two q — states with long range displacement and 
spin order dominated by E ff or E u phonons 3 . For both 
phonons there are 4 short bonds of length 1 — J'/k, and 
two long bonds of length yl+ 2J'/k + 5(J'/fc) 2 on each 
tetrahedra. The difference between the phonons is the 
arrangement of the short and long bonds on the lattice. 
In our simulations we use the q = state with E u phonon 
where all tetrahedra of the same orientation on the lattice 
distort in the same way. 

The simulations are performed on a pyrochlore lat- 
tice of 8788 spins. In order to allow for lattice defor- 
mation, yet maintain the essence of the spin frustration, 
the boundary conditions are periodic for the spins and 
open for the coordinates^. The energy minimization 
was done using a zero temperature Metropolis algorithm 
which stopped when the change in the simulated system's 
energy was below the machine precision of eight digits. 
The value of J' is variable while J = 1 and k = 10 for 
all runs. Further details concerning the simulation can 
be found in Ref^i. The error bars are determined from 
an evaluation of the finite size effects: the energy of the 
q = state was calculated using the simulation for differ- 
ent lattice sizes, and seen to converge to the theoretical 
value as the lattice size goes to infinity. Multiple runs 
of the simulation with the same initial conditions showed 
other computational errors to be negligible. It should be 
noted that despite the fact that Eq. \E\ is valid only for 
small fry, in the simulations fry can assume any num- 
ber. 

In Fig. |3 we present the energy per spin as obtained 
for 4 different calculations/simulations: (1) and (2) The 
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FIG. 3: The minimum energy per spin obtained by the simu- 
lation program as described in the text. In the inset we depict 
a tetrahedral unit of the pyrochlore lattice with its spins con- 
nected by springs to enable visualization of the magnetoelastic 
term. 



energy of the q — state given by the computer program 
before and after energy minimization respectively. (3) 
The final energy of the system after energy minimization 
of a disordered structure; neither the distortions nor spins 
entered an ordered state in this case. (4) Eq. adapted 
to a finite lattice. 

At low values of J' 2 / ( Jk) all four cases agree with each 
other as expected. However, for high J' 2 /(Jk) the simu- 
lation depended strongly on the initial simulation state, 
namely, the ground state manifold for the Hamiltonian is 
non-ergodic and upon cooling the system can enter one of 
many metastable states. The discrepancy between cases 
(1) and (2) stems from large values of fry in the simula- 
tion that are not consistent with the theory. To summa- 
rize the simulations: (I) the computer is unable to find a 
better state than the q = one; (II) the theory is valid 
in the range J' 2 /k < 0.1J; and (III) in this range the 
disordered states are very close in energy to the q = 
state. This proximity in energy between ordered and dis- 
ordered states means that low temperatures are needed 
to isolate the q = state. Eq.|S| suggests that this energy 
scale is set by J' 2 jk. 

Next we examine the temperature dependence of the 
correlations by simulating a slow reverse annealing pro- 
cess starting from the q = state. In order to follow the 
displacement and spin correlations we performed com- 
puterized magnetic and nuclear neutron scattering, ac- 
cording to Ref^, in the 111 direction. The numerical 
scatterings of several simulation states are shown in the 
inset of Fig. 0] In panels (a) and (b) we present a ref- 
erence scattering at zero temperature from the perfect 
pyrochlore lattice with Heisenberg spins and no magneto- 
elastic term. In panel (a) we see two peaks which are the 
trademark of the perfect pyrochlore lattice where the dis- 
tance between two planes of the same kind (kagome or 
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FIG. 4: Nuclear and magnetic computer-simulated neutron 
scattering in the 111 direction. Inset: (a) Nuclear scatter- 
ing from a perfect pyrochlore lattice, (b) Magnetic scattering 
from one of the ground states of the Heisenberg Hamiltonian 
(without magneto-elastic term) on this lattice, (c)-(f) Nu- 
clear and magnetic scattering from the states obtained by a 
computer simulation of the Heisenberg Hamiltonian with the 
magneto-elastic term with k = 10 and J' = 0.1 (Eq. |SJ at 
different temperatures (note the scale difference). In (c) and 
(d) T = 10~ 6 J and in (e) and (f) T = 10~ 4 J. Main figure: 
magnetic scattering amplitude as a function of temprature. 



triangular) is 2^/2/3i? and between two different planes 
(kagome-triangular) is ^/2/3-R. In panel (b) we see the 
lack of any spin correlations as expected from the Heisen- 
berg ground states. In panels (c) and (d) we show a 
similar experiment, but this time at finite temperature 
(T = 10~ 6 J) and with a magneto-elastic term. Un- 
der these conditions long range magnetic order is found. 
However, no effect is detected in the nuclear scattering 
within our resolution. Panels (e) and (f) are the same as 
(c) and (d) but at a higher temperature. This reduces 
both the nuclear and magnetic scattering due to ther- 
mal vibrations, but it affects the magnetic scattering to 
a greater extent. This magnetic scattering amplitude is 
plotted in Fig.0|as a function of temperature. Clearly it 
decreases as the temperature increases and is not mea- 
surable at T > J' 2 /k. This is the main result of the 
simulation. 

We use our numerical findings to determine the Hamil- 



tonian parameters of YMO in the framework of the 
magneto-elastic coupling. We compare the theoretical 
change in bond length, which, to first order, is set by 
J'/(kR) to the 5% distribution in Mo-Mo distance ob- 
served by Booth et al&. We also estimate J' 2 /k by the 
temperature at which spin freezing is observed in YMO, 
namely, 22.5 K. From these numbers we conclude that 
J' ~ 0.01 eV/A and k ~ 0.1 eV/A 2 . These numbers 
should be considered only as order of magnitudes; if one 
would take the melting temperature as the point at which 
the magnetic amplitude drops to 1/e of its T = value, 
for example, the estimate of both k and J' would increase 
by a factor of 10. For comparison in ZnCi'204, another 
well studied frustrated Heisenberg spin system with ap- 
proximately the same spin-spin bond length, J' = 0.04 
eV/Ai£ and k = 6.5 eV/A 27 . However, the reason YMO 
maintains its cubic symmetry at all temperatures is still 
a puzzle. One possibility emerging from the simulation is 
the lack of ergodicity; upon cooling the system is stuck in 
a metastable state and does not relax to its ground state. 
In this case the bulk cubic symmetry is maintained, but 
locally the system is under immense strain and slightly 
distorts from its ideal cubic structure. 

To summarize, the width of the internal field distribu- 
tion A, as detected by the /iSR, grows upon cooling at 
a rate which cannot be explained by the increasing sus- 
ceptibility alone. Therefore, we conclude that the width 
of the distribution of coupling constants 5 A also grows 
upon cooling. We attribute this to temperature depen- 
dent lattice distortions. However diffraction techniques 
show no structural changes and only local probes have 
revealed a distribution of bond lengths. To address this 
discrepancy we performed numerical simulations of the 
Heisenberg Hamiltonian with magneto elastic coupling 
at finite temperatures. This allowed us to estimate J' 
and k in Eq. and to suggest that upon cooling the sys- 
tem is trapped in a metastable state. 
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